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1. Introduction 

Recently there has been a renewed interest in the finite temperature dynamics of the 
one dimensional spin 1/2 Heisenberg model, especially on the question of diffusive spin 
transport]!], |2], |3|, ||. In particular, it was argued that the integrability of the model 
implies pathological spin dynamics and presumably the absence of spin diffusion H, |J. 
The role of conservation laws was pointed out in reference |7[] were it was shown that in 
several quantum integrable models the uniform (q = 0) current correlations do not decay 
to zero at long times. This result, established using the Mazur inequality ||, suggests 
pathological finite temperature dynamics. 

As far as the Heisenberg model is concerned, the analysis of conservation laws has 
shown that the energy current operator commutes with the Hamiltonian, suggesting 
anomalous finite- (g, u) energy density correlations. However, for zero magnetic field, 
this method turned out to be insufficient for deciding about the decay of the uniform spin 
current correlations. This case is closely related to the behavior of the finite temperature 
conductivity in the one dimensional model of spinless fermions at half-filling interacting 
with a nearest neighbor interaction (the "t-V" model 0). 

In this work, we address the issues raised above by the numerical diagonalization of 
the Hamiltonian matrix on finite size lattices. More precisely, we study the implications 
of the energy current conservation on the (q, u) energy density correlations, and, as an 
alternative route to the analysis of spin diffusion, we investigate the decay of the uniform 
(q = 0) spin current correlations. 

The paper is organized as follows: in section 2, we recall the Heisenberg Hamiltonian 
and define the various quantities studied below. In section 3 we briefly summarize the 
phenomenological picture of diffusion. There, we also argue that the decay of the uniform 
spin current correlations to a finite value is incompatible with a diffusive behavior, 
assuming continuity in the wave- vector q of the correlations at q — 0. Next, we test these 
ideas in section 4 in the XY limit, where results can be obtained analytically. Turning 
to the numerical results, in section 5.1 we present the energy density correlations at 
infinite temperature for the case of the isotropic Heisenberg model. A simple ansatz 
for the observed behavior suggests a logarithmic dependence at low frequencies for the 
energy autocorrelation function. As far as the spin dynamics is concerned, numerous 
studies of the (q,w) spin density correlations exist |], |], |3|. Therefore, in section 5.2, 
we restrict ourselves to the decay of the uniform spin current correlations for various 
values of the anisotropy parameter A and temperatures. Interestingly, it turns out that 
these do not decay to zero for A < 1. According to the argument given in section 3, 
this result implies non-diffusive spin transport. Section 6 contains a short discussion on 
experimental relevance of these findings and open questions. 
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The model 



The anisotropic Heisenberg Hamiltonian for a chain of L sites with periodic boundary 
conditions is given by: 

L L 

# = X> = 3 E^+i + SfSf +1 + AS?S? +1 ), (1) 
i=i i=i 

where Sf = \of-, erf are the Pauli spin operators with components a = x, y, z at site I. 

For a conserved quantity A = Y^t=i a ii [A,H] = 0, the continuity equation in 
q— space defines the current j q : 



^ = 2isin( g /2)j, (2) 



with 



a q 



= TfE ei ^> U=^fY,<2*3i (3) 

and a q (t) = e iHt a q e~ iHt . 

Setting ai = Sf, hi we find the following spin and energy currents respectively: 

i\ = J(S?S? +1 - S?S? +1 ) (4) 

•H T^(Q X QZQV CV qzqx \ 

Jl — J \ D l-l D l °l+l °l-l D l 

-t- J — O i _ 1 d l O l+1 ) 

+ J'AiS^SfS^ - Sl.SfS^) (5) 

For the discussion of dynamic correlations at finite temperatures, we chose to analyze 
the anticommutator form: 

S A A(q,t-t')= 1 -({a q (t),a_ q (t')}) (6) 

where ( ) is the thermal average at temperature T = 1/(3 over a complete set of states. 
Further, the frequency dependent correlation function defined by: 

/+oo 
duje^S AA (q,t) (7) 
■oo 

is symmetric in frequency, Saa(q, w) = Saa(q, —oS). 

A central point in our approach is the relation between the dynamic correlations of 
a quantity A and its corresponding current correlations, which we obtain by using the 
continuity equation @: 

uj 2 S AA (q,uj) = Asm 2 (q/2)S j A j A(q,u) (8) 



4 



In particular, we will discuss the asymptotic value of the current correlations 

<W = Km / jA r (9 = M) (9) 
J 3 t-*oo S jAjA (q = 0,t = 0) V ; 

A finite value of Cjaja translates to a 5(u) peak in SjAjA(q = 0,u) and, as we will 
discuss below, implies restrictions in the behavior of Saa(q, 

An important observation is that the energy current j H of the Heisenberg model 
commutes with the Hamiltonian|7j, so that CjHjH = 1, whereas the spin current does 
not. However, it will turn out that Cjzjz > for A < 1, meaning that the spin 
current and energy current correlations are similar in the sense that in their frequency 
representation, they both exhibit a finite weight S(u) function. 



3. Diffusive behavior 

When we consider the (q, u;)-dependent correlations of a conserved quantity A such as 
the magnetization, it is usually assumed, largely on phenomenological grounds, that 
they exhibit a diffusive behavior in the long-time \t — 1'\ — ► oo, short wavelength q — > 
regime 0: 

Wg,t-t')~e _iWM ' 1 (10) 
where Da is the corresponding diffusion constant, or 

2D A q 2 

for u) — > 0. 

This Lorentzian form correctly reduces to a 8{uS) function in the limit q — > 0, as 
implied by [A, H] = 0. Further, using the continuity equation (|]) for q — > 0, we obtain: 



2D A u 2 



SjA jA (q,u) ^ /n ^ ; 5 (12) 



which gives the diffusion constant .D^ when first, the limit q — > and then, cj — > are 
taken. On the other hand, if the current correlations for q = do not decay to zero 
at long times, CjAja > and SjAjA(q,u) has a finite weight 6(u) component which is 
incompatible with the diffusive form (0). In this reasoning, we must assume a regular 
behavior of the correlation functions in the q variable. 

To summarize the argument, if a quantity A is conserved ([A, H] = 0) and its current 
j A is either conserved {\j A , H] = 0), or Cjaja > 0, then continuity in q at q = excludes 
a diffusive form ( |TUD for the corresponding correlation Saa(<1, t — t'). 
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4. XY limit 



A simple model for testing these ideas is the XY limit (A = 0), of the Heisenberg 
model. In this case, both the energy current j H and the spin current j z commute with 
the Hamiltonian. The model can be mapped to a free spinless fermion model by using a 
Jordan- Wigner transformation which allows us also to evaluate explicitly the spin and 



energy dynamic correlations at /3 — 0. In the spin case, these are well known results flO 



ShhM = '"'t 12 '!:^' 72 0(\2J si„(f)| - M) (14) 
8 sin (|) 2 

These forms are indeed consistent with the conservation of both spin (energy) and spin 
current (energy current) as they reduce to a 5(u) function when the limits q — > 0, uj — > 
are taken. 

Further, the time decay of the autocorrelations at (3 = is not of the form 1/ y/i, as 
predicted by the diffusion hypothesis. Indeed, 

(smsn = \j 2 (Jt) as) 

(W = J {Jl{Jt) + Jl{Jt)) (16) 
which both behave as 1/t for t —>■ oo. 



5. Anisotropic Heisenberg model 

5.1. Energy correlations 

As we mentioned earlier, the energy current j H associated with the anisotropic 
Heisenberg model ([!]) commutes with the Hamiltonian for all values of the parameter 
A. Therefore, the time correlations do not decay at all (CjHjH = 1) and according to 
the argument explained in section 3, no diffusive energy transport occurs. However, the 
conservation of j H does not provide us with any details about the shape of Shh(o,uj) 
at finite q. In the absence of an analytical solution, we investigate this quantity by 
numerical diagonalization of the Hamiltonian matrix on a ring of 16 sites. 

In figure 1, we show Shh{<1i for A = 1, which is experimenally the most interesting 
point as it describes isotropic quasi one-dimensional antiferromagnets. We study the 
high temperature limit (3 = 0, which is the most convenient for a numerical study as it 
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involves the full excitation spectrum, but is also relevant experimentally for spin systems 
as the magnitude of J can be of the order of T. The plot is represented as histograms 
of width 0.06u;/J, all the frequencies which fall into one interval are summed up. The 
inset shows the normalized, integrated (prior to summing nearby frequencies) quantity 

fo + du/SHH(q,u') 

I H H{q,uj) = oo (17) 

Jo+ du'S HH {q,uj') 

which has the advantage of smoothing out the finite size discontinuities. To point out 
the practically linear integrated behavior of the pure Heisenberg model, we also show 
the same quantity for a more generic case obtained by adding a next-next neighbor 
(non-integrable) interaction J 2 - 

The simplest way to describe this behavior is by means of "plateaus" given by the 
following ansatz: 



Shh ^ w ) = V* I - <V3(1 " cosfa)) (18) 

loy 1 — cos(g) 



which satisfy the first J du Shh(q,w) = 3nJ 2 /8 and the second / du u 2 SHH(q, = 
37rJ 4 (l — cos(g))/8 exact moments for (3 = 0. Further, this ansatz is compatible with 
the limit Snniq — ► 0,u;) — > S(u) as implied by the conservation of energy. Using the 
continuity equation (|]), we obtain for small q and u 

s^Aw) = y ^ L ^ e ^\- ]ff J \i\) ( 19 ) 

which correctly reduces to a 5(a;)-function for q — > 0, in agreement with the conservation 
of the energy current [j H , H] = 0. 

Using this ansatz we find for the energy autocorrelation function (obtained by 
integration over q): 

(hi(t)h l )e iujt dt = C -C 1 \n(Lu/J) + O(LU 2 ) , C ,Ci>0 (20) 

a logarithmic behavior at low frequencies, in contrast to the diffusion form 

We should stress that these results are only indicative, as they are obtained from 
small size lattices which can provide reliable information only for correspondingly high 
frequencies and large wave-vectors. Nevertheless, the consistency of these results with 
the arguments presented above against a diffusion form are encouraging. 

5.2. Spin correlations 

The spin density dynamic correlations Ss^s z (q^) have been the subject of many studies 
which have not been able to answer the question of spin diffusion unambiguously. Here, 
we revisit this problem by investigating the compatibility between spin density and 
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spin current correlations, which requires that we calculate Cj*j*. In contrast to the 
energy current, the spin current j z does not commute with the Hamiltonian, so that 
Sjzjz = 0, u) is different from a pure S(u) function. Nevertheless, if Cj*j* > 0, which 
means that Sjzjz(q = 0,u) has a finite weight 5— function at uo = 0, our previous 
arguments against diffusion still hold. 

In determining Cj»j* , we noticed a peculiar difference in the low frequency behavior 
of Sjzjz(q = 0,u) depending on the anisotropy parameter A. In figure 2, we show 



the corresponding integrated, normalized quantity. We see that for A = cos(7r/n), 
n — 3, 4, . . . {n — 3 in the figure) all the low frequency weight of Sjzjz[q = 0,uj) is 
concentrated in the 5-function at u — 0. In contrast, for neighboring values such as 
A = 0.45 or 0.55, we observe a shift of weight to a low frequency region whose size 
decreases as the system grows (inset) and eventually vanishes as L — > oo. We believe 
that the behavior of this special A points is related to the existence of finite length 
strings (bound states) as they appear in the formulation of the thermodynamics of 
the Heisenberg model, within the Bethe ansatz method ^] . It seems that in order to 
determine Cjzjz from finite size systems for A ^ cos(ir/n), we should include the weight 
from these low frequency regions. As an example, doing so for A = 0.45 gives us a 
value of Cjzjz = 0.66 for L=16 (figure 2). Having discussed this technical issue, we can 
then determine Cjzjz for different size systems, as a function of temperature and A. 
By extrapolating our finite size results to the thermodynamic limit using second order 
polynomials in 1/L for L = 8, • • • , 18, we obtain the results shown in figure 3. Their 
striking feature is that for T > J, Cjzjz is finite in the A < 1 region, and practically 
zero when A > 1. In this regime, according to our previous argument, we expect a non- 
diffusive behavior. 

Deciding about the behavior of Cjzjz for A > 1 at finite temperatures is rather 
subtle. The reason is that in the Heisenberg model, A = 1 corresponds to a point 
of change of symmetry, from easy plane to easy axis, accompanied by the opening 
of a gap. In the fermionic version of the model, the "t-V" model, it corresponds to 
a metal-insulator Mott-Hubbard type transition, with the charge stiffness changing 
discontinuously|T2|1 at zero temperature. We should note that this discontinuity is 
difficult to reproduce by numerical simulations on small finite size lattices, as the 
transition corresponds to the divergence of the localization length. Considering that at 
high temperature, Cjzjz behaves similarly to the charge stiffness in the "t-V" model [0] 
we understand why it is difficult to decide whether Cjzjz is greater than zero in the 
region A > 1 and T < oo. For the same reason, we cannot exclude that Cj*j* behaves 
discontinuously at A = 1. Nevertheless, it seems unambiguous that Cj»j* ~ for 
A > 1.5 and T > J. 




(21) 
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6. Discussion 



The results presented are of interest in recent experimental studies [13] of spin dynamics 
in quasi-one dimensional materials such as CuGeO^ and SrzCuOs. Particular attention 
should be paid to the unusually high value of the diffusion constant found in NMR 
experiments on Sr^CuO^^,, perhaps related to the integrability of the Heisenberg 



model as discussed above. Furthermore, our results on the behavior of energy density 
correlations are of interest in the interpretation of the quasi-elastic Raman scattering, 
related to magnetic energy fluctuations] 13, |T5| . We should emphasize that no diffusion 
form should be expected for the energy density correlations in the isotropic Heisenberg 
model with only nearest neighbor interaction. An eventual diffusive behavior should 
be attributed to next-nearest neighbor coupling, interaction with phonons or deviations 
from one dimensionality. Finally, the main unresolved issue in this work is a better 
understanding of the finite temperature spin dynamics at the isotropic point. 
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Figure captions 



Figure 1. Energy density correlation function Snniq,^) at /? = for A = 1, 
q = (27r/16)n, n = 1,...,4. The inset shows the normalized, integrated quantity 
I HH (q,u) for A = 1, J 2 = and J 2 = 0.2 J. 



Figure 2. Integrated q = energy current correlations Ijzjz(u)) for TV = 16, /? = 
and A = 0.45, 0.5 and 0.55. The inset displays Ij* j* (oj) for L = 12, 14, 16, 18 at inverse 
temperature (3 = and A = 0.2. 



Figure 3. Values of Cp-p as function of A and f3 obtained by extrapolating second 
order polynomials in l/L from results on systems of sizes L = 8, . . . , 18. Calculations 
were done for A = cos(7r/n), n = 3, 4, 5, 6, 7, 10 and A = 1.0, 1.1, 1.5. 
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